Comment on "Direct Mapping of the Finite Tem- 
perature Phase Diagram of Strongly Correlated 
Quantum Models" 

In their Letter Zhou, Kato, Kawashima, and 
Trivedi claim that finite-temperature critical points of 
strongly correlated quantum models emulated by opti- 
cal lattice experiments can generically be deduced from 
kinks in the derivative of the density profile of atoms 
in the trap with respect to the external potential, k — 
—dn(r)/dV(r). In this comment we demonstrate that 
the authors failed to achieve their goal: to show that 
under realistic experimental conditions critical densities 
n c (T, U) can be extracted from density profiles with con- 
trollable accuracy. 

When illustrating their proposal with numerical data 
(in Figs. 4 and 5), the authors take it for granted that (i) 
sharp features in n only come from the critical behavior 
(Fig. 4), and (ii) critical behavior in a trapped system 
necessarily results in a cusp in k (Fig. 5). Both assump- 
tions are wrong, and this invalidates the first-principles 
component of the work. 

(i) Sharp features in k do not necessarily originate from 
critical fluctuations. A relevant counterexample would be 
supcrfluid helium in a Dewar with the density gradient 
sharply peaked at the wall. The only direct simulation 
performed by the authors and presented in Fig. 4 falls 
in this category. Both sharp features correspond to ex- 
tremely large gradients of particles (near the trap perime- 
ter) and holes (near the Mott insulator phase) , when the 
particle/hole densities change by ~ 100% over one lattice 
spacing in the radial direction! In this case it is funda- 
mentally impossible to extract the critical concentration 
of particles, n c (T, U), or holes, 1 — n c (T, U), with control- 
lable accuracy since the critical region is simply absent. 

(ii) The four illustrations presented in Fig. 5 are based 
on the local density approximation (LDA) reconstruc- 
tion, not the results of simulations in the trap, and are 
thus misleading. In all cases the density changes by 
~ 100% at a distance of just a few lattice spacings, mean- 
ing that the size of the critical region is limited, and 
features in k must be dramatically rounded. The effect 
of finite-size rounding of critical singularities is routinely 
seen in Monte Carlo simulations but the authors do not 
address this problem at all. By the very nature of second- 
order phase transitions LDA inevitably fails because of 
a divergent correlation radius in the vicinity of a critical 
point, and this is precisely why the sharp features get 
rounded. The reader should not be misled by the sharp 
features in k because pronounced kinks in Fig. 5 are hy- 



pothetical and do not represent the actual behavior of k 
in a trapped system! 

To show how severe the rounding of critical singulari- 
ties in a trapped system is, and how it makes the precise 
determination of critical parameters from k nearly im- 
possible, we performed simulations of the system speci- 
fied in Fig. 5(B). The first-principles data are presented 
in Fig. [[J 

In sharp contrast with the LDA plots of Ref. [l| , we ob- 
serve no features for the published parameters that would 
allow one to define errorbars through the full width at 
half maximum (or any characteristic interval) of the ex- 
perimental curve. This brings us to the above-formulated 
conclusion. It leaves the question open whether it is fea- 
sible to extract critical parameters from the density 




FIG. 1: Density profile and its derivative with respect to the 
external potential (inset) for parameters of Fig. 5b in Ref. [lj. 
Our unit is the hopping amplitude t. Errors are on the order 
of the point size, or smaller. 
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